Nucleation and cap formation on symmetric metal nanocatalysts: A first step towards 
chirality-controlled single-walled carbon nanotube growth 
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Symmetric minima of surface potential energy of a nanocatalyst act as nucleation sites for chirally 
selective initial growth of single walled carbon tubes at low temperatures. The nucleation sites are 
sites of maximum coordination number of the adsorbed carbon. We show this using the five fold 
symmetry of a pentagonal pyramid of an icosahedron. Initial zigzag structure from nucleation sites 
results in formation of hexagons and pentagons that result in anomalous cap formation. Possible 
cap lift off mechanism is discussed. 

PACS numbers: 61.48. De, 81.10.Aj, 81. 16. He, 31.15.A- 



The controlled growth of single-walled carbon nan- 
otubes (SWCNTs), i.e. with predetermined diameter and 
chirality angle, remains one of the holy grails of nanotube 
research despite two decades of substantial progress in 
their synthesis QUI]. The chirality indices (n,m) char- 
acterize the diameter and the chirality angle, and deter- 
mine whether a SWCNT is a metal or a semiconductor. 
At present SWCNTs with a mixture of chirality indices 
are synthesized in a variety of ways [5] , and growing them 
with specific chirality as desired is crucial for their appli- 
cations in future-generation electronics [6j. At the core 
of the problem is the complex nature of the growth con- 
ditions. While it has been possible to observe some as- 
pects of the growth in-situ H] , the atomistic dynam- 
ics that takes place is beyond the reach of experimental 
observation. Computationally the non-equilibrium pro- 
cesses have been studied with different techniques includ- 
ing density functional, tight binding, and classical molec- 
ular dynamics calculations [91114]. These studies have 
shed light on the growth mechanisms in general, but few 
have looked into how the chirality and diameter are set 
during the growth. Reich et al showed in the limit of large 
catalysts how certain pre-assembled carbon caps of root 
growth are preferred during nucleation due to their epi- 
taxial relationship with solid catalyst surface [T5] . Zhu ct 
al, however, questioned the rigid epitaxial relation used 
in the model at the interface between the cap and the 
catalyst [IB]. The work of Gdmez-Gualdrdn et al also 
shows how the growing nanotube affects the shape of 
the catalyst [17] . Nevertheless the formation of the cap, 
and the size of the catalyst are agreed upon generally 
to be major deciding factors in chirality selection dur- 
ing the initial stage of growth. The exact way how the 
cap itself is formed, however, has not been shown and 
has been rarely attempted jTTJ, [TH] . In the commonly in- 
voked vapor-liquid-solid (VLS) model to show the initial 
stages of SWCNT growth [T5], the nanocatalyst is first 
saturated with carbon atoms and then the carbon atoms 
form a graphitic sheet on the surface and subsequently lift 
off with incorporation of pentagons. During the course 
of the formation of the cap and during the subsequent 
growth of SWCNT by the incorporation of more carbon 
atoms at the interface of the nanotube and the nanocat- 



alyst, it is argued that ring transformations occur that 
anneal defects into hexagons over the span of the growth 
time (2D] |5I] • This is a double aged sword in that while 
the annealing of the defects to form hexagons is desirable 
on one hand, on the other it shows that at such high tem- 
peratures of growth that exceed the Stone- Wales trans- 
formations, it will be hard to retain the same chirality 
and local variations are bound to occur. Recently, it has 
been suggested because of the chirality instabilities that 
occur during the growth that it may not be possible to 
achieve chirally controlled growth with nanocatalysts [22] 
and that an alternative route from nanorings without cat- 
alysts may be more realistic. In contrast, Sankaran et al 
claim increasing the percentage of semiconducting SWC- 
NTs to ninety percent in the yield of plasma enhanced 
catalytic chemical vapor deposition (PECVD) by tuning 
the composition of bimetallic catalysts [3] . 

In this letter we show the possibility of growing SWC- 
NTs with predetermined chirality from symmetric solid 
nanocatalyst at low temperature. To our knowledge, the 
symmetry of the nanocatalyst and its effect on chirality 
selection has not been addressed. We assume the nucle- 
ation sites are symmetric and initate the growth of sym- 
metric SWCNT with pre-determined chirality. To test 
this idea, we chose to work with very small sized nanocat- 
alyst (diameter < 1 nm) because such nanocatatlysts 
lead to diameter-controlled SWCNT growth [53] [23]. Re- 
cently it has been shown that carbon atoms grow from 
nanocatalyst into a nanotube rather than encapsulate it 
in a fullerene cage due to the enhanced energy cost of 
encapsulation the smaller the nanocatalyst size is [55] , 

We follow a different simulation approach since we as- 
sume multiple adsorptions taking place more or less at 
the same time on the reactive sites of the symmetric faces. 
Therefore the number of atoms of carbon nucleating ini- 
tially will be at least equal to the reactive sites. This 
is markedly different from the molecular dynamics ap- 
proach that is often used where one carbon atom or a 
dimer or a trimer is fed to the growth process in peri- 
odic intervals or the Monte Carlo scanning and selecting 
of one site at a time based on some criteria. Such se- 
quential simulations of the reactions lead to meandering 
chains of carbon leading to tubular formations of mixed 



rings of carbon [51 [TJJ H3] ■ If the nanocatalyst surfaces 
are exposed to carbon atoms at the same time, the se- 
quential simulations can result in different outcomes than 
reactions taking place simultaneously. 

The Vienna Ab-initio Molecular Dynamics Package 
(AMD) is used to perform the calculations with the 
projector augmented wave (PAW) potentials 26J. The 
AMD is computationally very intensive, so we chose Fei3 
as our nanocatalyst. Fei3 could be dismissed as unre- 
alistic small system, but it has been used to simulate 
the growth of SWCNT [TT], the adhesion strength of 
monoatomic catalysts of the iron family [30] and recently 
to study bimetallic catalysts [31]. From symmetry point 
of view of our work, it has similar properties as the larger 
more realistic Fes5 cluster and thereby one could study 
the underlying physics involved in icosahederal nanocat- 
alysts without introducing possible complexities associ- 
ated with larger clusters [ST] , 

The AMD is executed on a picosecond time range while 
the experimental time is on the order of nanoseconds. 
Therefore, it is criticized as a method that cannot cap- 
ture the experimental processes that yield the final out- 
come. In this work we argue that the rings and the chiral- 
ities are set at the initial stages of growth, particularly 
in the low temperature range of 100 - 500°C. Cantaro 
et al have grown SWCNTs with thermal chemical vapor 
deposition (CVD) at temperatures as low as 350°C [28] . 
The possibility of lower temperatures and activation en- 
ergies than thoses of thermal CVD by PECVD has been 
suggested [251 . 

We performed ground state calculations first [57] . The 
resulting optimized structures were then heated to finite 
temperatures by scaling the velocities with temperature. 
Between the temperature steps, a micro-canonical ensem- 
ble is simulated. This approach is much faster than the 
computationally intensive simulation of canonical ensem- 
ble. 

The icosahedra structure of Fei3 is highly symmetric. 
Each atom of Fei3 forms one of the apexes of a pentag- 
onal pyramid when observed from the top of it. There 
are three possible adsorption sites for carbon atoms: the 
apex, the bridge between the apex and one of the atoms 
at the base of the pentagonal pyramid and the triangu- 
lar face of the pyramid. The adsorption energies were 
computed for the three possible sites at T = K first. 
The carbon placed on the bridge was pushed to the side 
and got absorbed on the face and resulted in the lowest 
energy. The triangular face allows the carbon atom to 
bond with the maximum number of Fe atoms (three in 
this case) and thereby results in lowest energy. The ad- 
sorption energy is equal to 7.39 ev which is in good agree- 
ment with that of [llj . The relative adsorption energies 
for the sites with respect to this energy are given in Ta- 
ble I. The question, however, is that are these results for 
ground state applicable for the finite temperature growth 
of SWCNT? To test their validity, we heated Fe i3 to 400 
- 800 K and found out the icosahedral structure of Fei3 
remained unchanged. Therefore, the carbon atom will 



see the same surface potential. To check whether the ki- 
netic energy of the carbon is large enough to overcome 
the local minima or the diffusion energy barrier on the 
surface when it moves into the potential well, we heated 
the optimized system, Fei3 + C, to 400 - 800 K as well 
and the results are given in Table I. The adsorbed car- 
bon on the face of the optimized structure remained on 
the face for both cases from T = K, with energy differ- 
ence of 0.14 to 0.52 ev between the two structures. On 
the other hand, the carbon attached to the apex became 
slanted at 400 K and 600 K and finally ended over the 
center of the face at 800 K. These results showing that 
the carbon is still bonded to the face at low tempera- 
tures, despite its kinetic energy, validate our assumption 
that the adsorption sites of T = K still hold for low 
temperatures. 

In our model, carbon atoms flow onto the face from 
one direction, for instance in PECVD. After placing the 
carbon atoms on the symmetric faces, we determined the 
next likely adsorption site by placing carbon atoms in 
three different places: i) along the bridge again (the car- 
bon atom bonds with four atoms, two from the vertices 
of the triangular face and two from the adsorbed carbon 
atoms on adjacent triangles), ii) on the apex again, and 
iii) on top of one of the adsorbed carbon atoms on the 
face. During the adsorption, the carbon atoms tend to 
redistribute themselves or spread out as much as possi- 
ble to reduce the strain energy. Therefore to promote 
unidirectional growth of nanotubes a geometrical con- 
straint has to be imposed to prevent the spreading out 
of the carbon atoms in all directions on the catalyst sur- 
face. Such an assumption has been applied previously 
by passivating the lower portion of the nanocatalyst [llj . 
We rule out the spreading out of carbon atoms on Fei3 
icosahedral surface into region below the Fe atoms form- 
ing the pentagonal base. The optimized structures of the 
second and third cases violate the geometrical constraint 
and were ruled out. The first case resulted in the lowest 
energy. The carbon atom is pushed outward along with 
the adjacent carbon atoms adsorped on triangular faces, 
as seen in the Fig. 1(a). This outward push, however, is 
counterbalanced by adsorption at diametrically opposite 
site when multiple adsorptions occur, and the push out- 
ward is not as significant as it looks in the figure. In the 
process, it forms the portion of the zigzag structure that 
will eventually form around the pentagonal base as seen 
in Fig. 1(b). This configuration is unstable energetically 



TABLE I. Relative adsorption energies (ev). Finite tempera- 
ture energies are obtained by heating the ground state struc- 
tures. Initial adsorption site is where carbon is placed before 
optimization. 

Initial Adsorption sites OK 400 K 600 K 800 K 



Face 
Bridge 
Apex 



0.07 
0.00 
2.13 



0.67 
0.82 
2.51 



1.09 1.57 

1.45 2.09 

2.52 1.14 



and if left alone to relax would optimize into a relaxed 
structure shown in Fig. 1(c). However, in the presence 
of carbon atoms in the vicinity, the zigzag structure with 
its kinks is more reactive and forms hexagonal rings as 
seen in Fig. 1(b). 
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FIG. 1. (color online) a) Formation of zigzag fragment, b) 
Initial hexagonal rings formation, c) Open structure from 
relaxed zigzag ring. The T^ K structures were obtained by 
heating the T = K ones. 

At T = K we supplied carbon in the form of 
uniformly distributed carbon atoms, carbon atoms + 
dimers, dimers only, and a ten-atom carbon ring. Three 
of the distributions, carbon atoms, carbon atoms + 
dimers, and the ring resulted in hexagonal ring forma- 
tion. The ten-atom ring when heated to finite tem- 
perature caused the distortion of Fei3, therefore, we 
proceeded only with atoms, and atoms + dimers. For 
brevity, the result for atoms only is shown in Fig. 1(b). 
One could raise the question of supplying randomly dis- 
tributed carbon atoms instead of uniformly distributed 
ones, but we argue in the tiny diameter range of Fei3 
there is not much volume for randomness and a symmet- 
ric supply is a good approximation. The relaxed struc- 
ture of Fei3 + IOC was also exposed to the same car- 
bon distribution at T = K, optimized and then heated. 



The optimized structure were open and are shown in Fig. 
1 (c) . The result underscores the fact that a zigzag struc- 
ture with kinks at the onset is critical to the formation 
of hexagonal rings and the icosahedral geometry is con- 
ducive to such a formation. The icosahedral structure of 
the Fei3 remained intact except for the local distortions 
in the vicinity of the adsorption sites at the interface. 
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FIG. 2. (color online) At T = 600 K a) 10-atom ring repulsed 
from Fei3 + C20 cluster, b) hemispherical C20 cap formation. 

In Fig. 2, we show the results of simulation of addi- 
tional carbon atoms. Results for atoms + dimers are 
similar therefore, are not shown for brevity. In both 
cases of the input, a ring as shown was formed and was 
pushed off instead of forming part of the existing struc- 
ture. There are two triangles formed at the edge of the 
Fei3 + 20C on opposite sides as neighboring dangling 
bonds combined. We broke the bond forming the trian- 
gle for each triangle and optimized again with the addi- 
tional ten carbon atoms on the top. The 10-atom ring 
was once again formed and was repulsed. At such low fi- 
nite temperature, the kinetic energy of the carbon atoms 
is not strong enough to overcome the repulsive barrier 
to attach to the existing structure. This could be ex- 
plained by the recent findings of charge transfer from the 
nanocatalyst to the carbon atoms |32J . The carbon atoms 
adjacent to the iron atoms become negatively charged, 
and this first set of negatively charged carbon atoms at- 
tracts the next batch of neutral carbon atoms by polar- 
izing them and integrates them into the structure. This 
polarization effect pushes the electrons to the edge of the 
nascent tube. These electrons at the edge repulse the 
7r electrons of the ring of Fig. 2 that was formed dur- 
ing the interaction. The repulsive force dominates over 
whatever remnant attractive force that exists from the 
dangling bond of the edge carbon atoms trying to bond 



with other carbon atoms on the edge [33J and the at- 
tractive force emanating from the iron atoms, resulting 
in the pushed away carbon ring in Fig. 2. We reduced 
the number of carbon input atoms to five instead of ten 
to reduce the repulsive force and to see if the attractive 
force could lead to a closed structure instead of an open- 
ended structure with dangling bonds. The two triangles 
at the edge were removed by breaking the bond between 
the neighboring carbons at the edge before placing the 
five carbon atoms [33]. This leads to the hemispherical 
cap of C20 fullerene shown in Fig. 2. The cap is not 
the usual hemispherical Cgo cap that is often attached 
to the end of the nanotubes, and does not obey the iso- 
lated pentagon rule. This anomaly has been recently 
shown to exist at the end of the (3,3) SWCNT [3? . The 
ultra-narrow SWCNTs require larger energy of curvature 
and therefore more pentagons that are not isolated are 
necessary to form a cap. This is counter to the conven- 
tional isolated pentagon rule, but may be feasible for the 
ultra-narrow SWCNTs. The larger energy of curvature 
for these ultra-narrow SWCNTs can lead to well defined 
chirality in contrast to larger diameter SWCNTs whose 
energy of formation are almost identical |15| leading to 
a mixture of chiralities. Moreover, in an entirely differ- 
ent context it has been shown that SWCNTs grown from 
nanodiamonds have polyhedron-shaped caps rather than 
the usual hemispherical ones [35 . Thus the formation of 
a Cgo type cap may not be the necessary condition for 
the narrowest SWCNTs. The forces that are available at 



the nucleation: the dangling bonds trying to form bonds 
with their counterparts and the strong adhesion energy 
of monometallic catalyst such as iron could provide the 
curvature energy necessary. On the other hand as noted 
above the repulsive force acts as a countervailing force 
against the attractive force and possibly leads to the cap 
lift-off. 

In summary we have shown nucleation on reactive sites 
of symmetric nanocatalyst could lead to a well defined 
chirality of ultra-small SWCNTs at low temperature with 
icosahedral Fei3. The small system results in well defined 
chirality because of the relatively large energy of curva- 
ture. A hemispherical cap of C20 is formed which could 
lead to the growth of (5,0) nanotube. An electrostatic re- 
pulsion counteracting with attractive forces of dangling 
bonds and the nanocatalyst is responsible for cap lift-off. 
Plasma enhanced chemical vapor deposition can supply 
the carbon atoms or carbon atoms + dimers which can 
be incorporated through root growth mechanism for ex- 
tended growth. 
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